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RANDOM  SEQUENTIAL  PACKING  IN  EUCLIDEAN  SPACES 
OF  DIMENSIONS  THREE  AND  FOUR 
AND  A  CONJECTURE  OF  PALASTI 

By 

B.  Edwin  Blaisdell  and  Herbert  Solomon 


INTRODUCTION. 

Random  sequential  packing  of  hyper cubes  in  a  larger  hypercube 
and  a  conjecture  of  Palasti  [11]  that  the  limiting  density  8^  in 
a  space  of  dimension  d  equals  6**  where  8  is  the  limiting  packing 
density  in  one  dimension,  i.e.,  on  a  line,  continues  to  be  studied, 
but  with  inconsistent  results.  Renyi  [12]  and  Dvoretzky  and  Robbins 
[5]  have  derived  an  integral-difference  equation  for  g  and  Blaisdell 
and  Solomon  [3],  in  an  investigation  of  the  Palasti  conjecture,  have 
calculated  the  value  of  B  to  14  significant  digits,  for  purposes  of 
this  paper  approximately  0.7476.  In  that  paper,  we  have  also  studied 
packing  on  finite  lattices  with  rigid  boundaries  by  computer  experiment 
in  one  and  two  dimensions.  In  one  dimension  our  results  were  in  good 
agreement  with  our  extension  of  an  analytic  result  of  Mackenzie  [10] 
for  lattice  lines,  namely 

(1)  c  -  1- tl  +  a/n  +  0(2n/2/(n/a).,nk0-k1a_1  +  0(a"2)] 

where  B  “  11m  c  as  an  ^  0  and  a  ^  0  » 

and  kQ  -  1-3,  \  -  ^[3-exp(-2Y)l  ~  0.2162  and  Y  is  Euler's 
constant,  n  Is  the  number  of  lattice  nodes  on  the  bounding  lattice 
line  and  a  Is  the  number  of  lattice  nodes  on  the  edge  of  the  randomly 
deposited  lattice  segment.  Define 
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(2) 


xx  *  (l-c)/(l+a/n)  -  0.2524  +  0.2162a 

by  a  partial  rearrangement  of  (1).  This  permits  satisfactory  extra¬ 
polation  of  our  finite  one-dimensional  lattice  data  to  the  infinite 
continuum.  A  least  squares  treatment  of  our  one-dimensional  data  in 
[3]  yielded 

(3)  xx  z  0.0000(+  0.0001)  +  0.0038 (+  0. 0036)a~2 

where  the  values  in  parenthesis  are  the  estimated  standard  deviations 
of  the  preceding  estimated  parameters.  Note  that  the  estimated  inter¬ 
cept  is  zero,  thus  agreeing  with  equation  (1)  (that  is,  the  lattice 
lines  situation).  A  similar  treatment  in  [3]  of  our  two-dimensional 
data  yielded 

(A)  x2  -  xd  «  -  0.0025 (+0.0001)  +  0.0109 (+0.0023) a-2 

where  is  the  same  as  in  equation  (2)  except  that  c  is  re¬ 

placed  by  (cd)1/d.  The  subscript  of  x  Indicates  the  dimensionality 
of  the  lattice  space.  This  result  Is  evidence  for  a  small  but  signi¬ 
ficant  discrepancy  when  Palastl's  conjecture  is  viewed  in  two  dimensions, 
l.e. ,  the  packing  density  is  larger  in  two  dimensions  than  Palastl's 
assertion. 

In  1975  Akeda  and  Horl  published  a  note  [1]  in  which  verification 
of  Palastl's  conjecture  In  the  plane  was  asserted  on  the  basis  of 
computer  experiments  in  two  dimensions  on  squares  with  n/a  ■  100  and 


a  -  the  value  of  the  floating  point  mantissa  of  their  computer ,  i.e., 
the  Continuous"  model.  After  being  informed  by  us  of  our  paper  [3] 
they  carried  out  further  experiments  [2],  performed  the  necessary 
extrapolation  to  a/n  -►  0  and  obtained  a  value  x^  *  -0. 0027 (+0. 0002)  . 
Akeda  and  Hori  also  studied  the  three-dimensional  "continuous"  model 
by  a  different  and  undescribed  method  for  values  of  a/n  •  8,  10*  12, 

15,  20  and  30  and  on  extrapolation  to  a/n  0  obtained  a  value 
x^  -  0. 0018 (+0. 0008) .  The  discrepancy,  x^,  has  a  sign  opposite  to 
that  for  two  dimensions  and  the  line  through  the  individual  points 
crosses  the  expected  extrapolation  line,  both  results  unappealing 
intuitively. 

In  1978  Weiner  published  a  purported  proof  of  the  Palasti  conjec¬ 
ture  [17]  but  the  validity  of  his  proof  has  been  challenged  by  several 
correspondents  [7,8,13,14,15, 16]  and  we  agree  that  his  claim  is  not 
valid.  Finegold  and  Donnell  [6]  published  a  note  in  1979  on  computer 
experiments  on  the  two-dimensional  "continuous"  model  with  periodic 
boundaries  by  a  "fine  mesh"  and  a  "coarse  mesh"  method  and  asserted 
verification  of  the  Palasti  conjecture.  Their  three  summary  values 
are  x£  «  -0. 0027 (+0. 0005) ,  -0. 0040 (+0.0005)  and  -0.0034 (+0.0012) . 

They  assert  that  the  use  of  periodic  boundaries  makes  it  unnecessary  to 
extrapolate  to  a/n  -►  0,  a  statement  which  is  only  true  in  one-dimension 
for  >  14.  The  accuracy  of  the  approximation  is  unknown  for  higher 
dimensions.  In  1980  Jodrey  and  Tory  [9]  published  the  results  of 
extensive  computer  experiments  on  the  "continuous"  model  with  periodic 
boundaries  for  dimensions  1,  2,  and  3  and  a  few  results  for  dimension 
4.  Their  extrapolated  values  are  *  0.0004 (+0.0005) , 
x^  ■  -0.0021  (+0.0002) ,  and  x^  ■  -0.0029(+0.0002)  and  the  average 
of  their  3  values  for  x^  -  -0.0003 (+0.0006) . 
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The  values  of  obtained  by  three  of  the  sets  of  authors 

[2,3,9]  are  In  good  agreement  -0.0027,  -0.0025,  and  -0.0021 
respectively  and  show  a  small  departure  from  the  Palastl  conjecture. 

The  values  of  obtained  by  two  of  the  authors  [2,9]  are  in 

poor  agreement  and  of  opposite  sign  +0.0018  and  -0.0029  respec¬ 
tively. 

Several  of  these  authors  have  also  reported  the  variance  of  the 
limiting  density  and  a  similar  account  may  be  given.  In  one  dimen¬ 
sion  our  results  were  in  good  agreement  with  an  analytic  result  of 
Mackenzie  [5]  for  the  lattice  line.  Define 

(5)  y1  -  s2/((a/n)(l+a/n) ]  +  A2(a) 

as  a/n  0  and  1/a  ■+  0  and  A2(a)  is  an  undetermined  power  series 
in  1/a.  A  least  squares  treatment  of  our  one-dlmenslonal  data  in  [3] 
yielded 

(6)  yx  -  0. 038K+0. 0003) +  0.016K+0. 0065)  (1/a) -0.0422(+0. 0238)  (1/a) 

where  0.0381  is  in  good  agreement  with  the  result  (0.038156)  we 
obtained  in  [3]  by  numerical  integration  of  the  analytic  expression 
for  the  continuous  line  obtained  by  Dvoretzky  and  Robbins  [5].  A 
plausible  generalization  to  d  dimensions  of  equation  (5) ,  namely 

(7)  yd  -  e^/[(a/n)(l+a/n)]d 
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TABLE  I 


Average  normalized  variances  for  sequential  random  packing  densities 

Ref 

Dimension 

yd 

notes 

(1) 

2 

0.0402  + 

0.0213  a 

(2) 

2 

0.0482  + 

0.0135  a 

3 

0.0413  + 

0.0206 

1 

0.0400  + 

0.0058 

(9) 

2 

0.0468  + 

0.0069 

3 

0.0451  + 

0.0101 

4 

0.0324  + 

0.0176 

(3) 

1 

0.0393  + 

0.0014  c 

2 

0.0526  + 

0.0027  c 

this  work 

3 

0.0519  + 

0.0062 

4 

0.0518  + 

0.0158 

Notes:  (a) 

only  the  values  for 

squares  of  edge  40  have  been  used. 

(c)  only  a  subset  of  the  values  obtained  In  (3)  have  been  used, 
as  described  below. 
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was  found  by  least  squares  to  be  in  good  agreement  with  our  two-dimen¬ 
sional  data  [3],  namely 


(8)  y2  -  0. 0508(+0.0010)+1. 3697 (+0. 3905) (l/a)3-3.6660(+l. 1400) (1/a)4  . 

Average  values  for  y^  are  given  in  Table  I .  The  factor  1+a/n 
which  normalizes  for  finite  hypercube  size  is  not  used  for  experiments 
with  periodic  boundaries. 

New  Developments 

We  have  now  conducted  experiments  in  three  and  four  dimensions. 

In  doing  so,  we  have  adhered  to  our  earlier  policy  of  studying  finite 
lattices  with  rigid  boundaries.  Although  this  policy  results  In  loss 
of  accuracy  when  extrapolating  to  the  infinite  continuum  [9],  it  is 
easier  to  program  and  permits  an  exact  accounting  of  every  lattice 
site  as  occupied  or  unoccupiable.  Authors  who  have  studied  the 
"continuous”  case,  l.e.,  have  used  floating  point  arithmetic,  have 
not  agreed  with  each  other  nor  with  our  present  results  in  dimension 
three.  One  possible  explanation  for  this  is  the  occurrence  of  holes  that 
present  a  very  tight  fit  in  at  least  one  dimension,  and  which  have  been 
missed  because  of  round-off  error  in  the  continuous  case.  The  likelihood 
of  these  tight  fits  occurring  will  increase  with  increasing  dimension  and 
may  account  for  Jodrey  and  Tory  [9]  finding  densities  increasingly  lower 
than  ours  in  dimensions  two,  three  and  four. 

Our  calculations  were  made  as  described  before  in  [3]  except  that 
we  now  used  a  DEC  PDP11/55  computer  running  FORTRAN  programs  under 


operating  system  RSX11M.  Lattice  points  were  identified  with  bits 
in  the  computer's  16-bit  words,  which  permitted  use  of  hypercubes 
containing  up  to  70  or  24  (about  330,000)  lattice  nodes.  The  new 
values  obtained  in  dimensions  three  and  four  are  given  in  Table  II. 

These  values  together  with  subsets  of  those  obtained  formerly  in  [3] 
in  dimensions  one  and  two  have  been  treated  together  since  a  regular 
pattern  related  to  dimension  has  become  apparent.  For  dimension  one 
we  used  values  for  a  £  100,  n/a  £  20  and  for  dimension  two,  values 
for  a  £  20,  n/a  £  10.  This  yielded  sets  of  about  the  same  number 
for  each  of  the  four  dimensions  so  that  the  large  number  of  values 
otherwise  available  for  dimensions  one  and  two  would  not  outweigh 
those  calculated  for  dimensions  three  and  four.  The  numbers  for  dimen- 

1 

sions  one,  two,  three,  and  four  are  respectively  28,  34,  51,  and  32 
including  all  values,  and  18,  22,  30,  and  24  excluding  values  for 
a  ■  2  and  n/a  35  5  which  we  found  as  before  [3]  to  require  more 
terms  in  the  expansion  of  equation  (1)  for  obtaining  a  satisfactory 
fit. 

We  have  examined  our  results  guided  by  the  theoretical  equation 
(2)  for  one  dimension,  substituting  for  c,  the  value  c^^,  (the 
Palasti  conjecture)  in  higher  dimensions.  The  results  of  our  examina¬ 
tion  are  given  in  Table  III.  Lines  1,  2,  3,  and  4  contain  the  least 
squares  results  treating  each  dimension  separately.  It  is  apparent 
that  the  intercept  (a  measure  of  the  departure  from  the  Palasti  conjec¬ 
ture)  increases  rather  regularly  with  dimension  and  that  the  coefficient 
2 

of  the  term  in  1/a  (a  measure  of  the  departure  from  the  limiting 
equation  (2))  decreases  with  dimension.  Lines  5,  6,  7,  and  8  contain 
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Lattice  Dimension  3  Lattice  Dimension 
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the  least  squares  results  treating  each  dimension  separately  but 
Including  least  squares  estimates  of  these  observed  trends.  The 
fit  is  of  course  the  same  but  the  near  constancy  of  the  respective 
parameter  estimates  shows  that  the  model  form  is  satisfactory. 

Lines  9,  10,  11,  and  12  contain  least  squares  results  treating  all 
dimensions  at  once  but  successively  adding  the  less  well  fitted 
values  for  a  *  2  and  n/a  *  5.  The  standard  error  of  estimate 
Increases  about  25%  but  the  F  value  increases  about  30%,  because 
of  the  substantial  increase  in  degrees  of  freedom,  as  these  values 
are  included.  A  search  was  made  for  further  terms  which  might  sig¬ 
nificantly  improve  the  fit  using  program  P9R,  all  possible  subsets 
regression,  of  the  BMDP  package  [4],  Lines  13,  14,  15,  and  16  con¬ 
tain  the  least  squares  results  for  the  best  choice  of  4  terms  from 
the  following:  (d-1),  1/a,  a/n,  1/a2,  1/a3,  1/a*,  l/(a2^d),  (a/n)(l/a2). 
The  additional  terms  make  a  significant  improvement  in  the  standard 
error  of  estimate  if  the  data  for  a  ■  2  or  n/a  *  5  are  included 
but  not  otherwise. 

The  residuals  from  the  fit  in  line  9  are  plotted  against  (d-1) 
in  Figure  1  and  against  l/(a2*/d)  in  Figure  2.  The  plotted  digits 
show  the  number  of  superposed  values  at  that  location.  There  appears 
to  be  no  trend  in  the  residuals  with  either  term,  further  testimony 
that  the  model  is  satisfactory. 

Computer  experiments  on  the  random  sequential  packing  of  finite 
lattices  with  rigid  boundaries  in  dimensions  one,  two,  three  and  four 
indicate  a  discrepancy  with  the  Palastl  conjecture  in  the  limit  l/a+0, 
a/n  +  0  which  Increases  about  0.0025  in  Per  dimension.  These 
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correspond  to  discrepancies  for  dimensions  two,  three  and  four  of 
0.0037,  0.0084,  and  0.0127  In  the  density  Itself,  c^  being  greater 
than  cd.  This  may  be  seen  by  a  rearrangement  of  the  limiting  form 
of  equation  (2)  with  c  replaced  by  c^d.  We  can  now  write 
cJ/d-Cl  *  z(d-l)  where  the  right  side  is  the  limiting  expression 
found  in  Table  III.  Then 

cd  -  [c1  +  z(d-l)]d 
or 

2  3 

cd-cj  -  zcj_1d(d-l)  +~  cd-2d(d-l)3  +  ^-  cd-3d(d-l)4(d-2)  +  ••• 

where  c^  -  .7476  and  z  *  i/c^-c^  ~  -0025,  the  excess  over  the  value 
given  by  the  Falastl  conjecture  when  d  *  2.  Note  that  the  expansion 
truncates  at  (d+1)  terms  and  because  z  -  .0023  is  quite  small,  no 
terms  beyond  the  first  are  significant  to  four  decimal  places  and  this 
leads  to  the  discrepancies  for  dimensions  two,  three,  and  four  given 
above. 
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